An Adaptive Algorithm for Ordinary, Stochastic and Partial Differential Equations

The theory of a posteriori error estimates suitable for adaptive refinement is well established. This work focuses on the fundamental, but less studied, issue of convergence rates of adaptive algorithms. In particular, this work describes a simple and general adaptive algorithm applied to ordinary, stochastic and partial differential equations with proven convergence rates. The presentation has three parts: The error approximations used to build error indicators for the adaptive algorithm are based on error expansions with computable leading order terms. It is explained how to measure optimal convergence rates for approximation of functionals of the solution, and why convergence of the error density is always useful and subtle in the case of stochastic and partial differential equations. The adaptive algorithm, performing successive mesh refinements, either reduces the maximal error indicator by a factor or stops with the error asymptotically bounded by the prescribed accuracy requirement. Furthermore, the algorithm stops using the optimal number of degrees of freedom, up to a problem independent factor.